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We propose an approach to enhance the mechanical effects of single photons in a two-mode 
optomechanical system. This is achieved by introducing a resonance-frequency modulation to the 
cavity fields. When the modulation frequency and amplitude satisfy certain conditions, the mechan¬ 
ical displacement induced by single photons could be larger than the quantum zero-point fluctuation 
of the oscillating resonator. This method can be used to create distinct mechanical superposition 
states. 
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I. INTRODUCTION 

One of the most interesting regimes in cavity quantum 
optomechanics [1, 2] occurs when the radiation pressure 
of a single photon is strong enough to push a mechanical 
oscillator with a displacement greater than the oscilla¬ 
tor’s zero-point fluctuations before escaping the cavity. 
In such a regime, one can generate distinct superposition 
states of a macroscopic object (the oscillator) [3, 4], and 
the quantum nonlinearity in the ultrastrong optomechan¬ 
ical coupling regime can lead to various quantum effects, 
including photon blockade [5-8], non-Gaussian states of 
the mechanical oscillator [9-13], modified spectra of pho¬ 
ton scattering [6, 14], and dressed-state-representation 
dissipations [15]. Specifically, this regime is defined by 
the condition: 

go > iOM > Kc, (1) 

where go is the single-photon optomechanical coupling 
strength, lom is the frequency of the mechanical oscilla¬ 
tor, and Kc is the cavity field damping rate. 

The requirement of go > ujm in (1) is understood from 
a generic single-mode time-independent optomechanical 
system with the Hamiltonian TL = uJcCi'a -I- ojMb^b — 
goa^a{b-\-b^), where a {b) and (6’!^) are the annihilation 
and creation operators of the cavity field (mechanical os¬ 
cillator). By solving the Schrddinger equation, one finds 
that the oscillator, when driven by a single photon, can 
evolve to a coherent state \2go/u}M) from the ground state 
after a time t = nfuiM- Therefore, go > wm provides a 
significant displacement away from the ground state. 

However, the realization of condition (1) has been chal¬ 
lenging for current experiments because go is small in 
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most optomechanical systems, and wm has to be suffi¬ 
ciently high to overcome thermal noise effects. In par¬ 
ticular, the quantum limit of phonon number nt, ob¬ 
tained by resolved sideband cooling is rib = ^c/(16w^) 
[16-20]. It should be noted that condition (1) is based 
on the single-mode time-independent Hamiltonian PI, 
and thus an interesting question is whether it is possi¬ 
ble to achieve ultrastrong coupling effects in multimode 
and time-dependent systems. Indeed, a similar ques¬ 
tion has been explored in electromechanical systems re¬ 
cently [21, 22], and it is found that a suitable modulation 
can be used to effectively enhance photonic nonlinearities 
in electromechanical systems. In addition, an interfero¬ 
metric scheme has been proposed to detect optomechan¬ 
ical coherent interaction at the single-photon level [23]. 

The methods of enhancing radiation pressure effects 
of a single photon in optomechanical systems have 
been recently discussed by several authors [24-28]. In 
this paper, we approach this problem by modulating 
the photon-tunneling process in the “membrane-in-the- 
middle” (MIM) configuration. The MIM optomechanical 
system has been studied experimentally [29-37] and the¬ 
oretically [38-44]. Here we show that by changing the 
frequency of the cavities in a time-dependent manner, 
the photon tunneling oscillation can be modulated. We 
discover a set of conditions such that a single photon can 
significantly displace the membrane, even though go I to m 
is much less than 1. Below we derive an effective Hamil¬ 
tonian of the system and present numerical evidence of 
the ultrastrong coupling effect. 


II. THE MODEL 

We consider an MIM optomechanical system (Fig. 1) 
in which a vibrational mode of the membrane is coupled 
to the left and right cavity field modes via the radiation 




FIG. 1. (Color online) Schematic of the optomechanical sys¬ 
tem with the “membrane-in-the-middle” configuration. The 
cavity helds couple with the mechanical oscillation via the 
radiation pressure coupling go, and the optical fields in the 
left and right cavities couple to each other via the photon¬ 
hopping interaction J. To realize the cavity-frequency mod¬ 
ulation in Eq. (3), the two outside mirrors could be forced 
to oscillate with the same displacements Xn(t) = Xu{t) = 
{Ao/oJc)L cos(2 Jt) by using an external device such as a piezo¬ 
electric transducer. 


pressure difference between the two cavities, and tunnel¬ 
ing of photons through the membrane is allowed. The 
Hamiltonian of the system [38-44] is modeled by {h = 1) 

H = WLO^az, -I- ujRal^an -|- J(a]^a_R -I- a^jai) 

+u!Mb'‘b +go^a'laL - al^aR){b + b'‘), ( 2 ) 

where and b are the annihilation operators of the 

left (right) cavity and the mechanical membrane, with 
the respective resonance frequencies uir(^r) and ujm- The 
parameters J and are the coupling strengths of the 
photon hopping and the optomechanical interaction, re¬ 
spectively. 

A key feature in our model is the introduction of a weak 
frequency modulation to the left and the right cavities 
such that wl and lor in Eq. (2) are given by 

WL(t) = Wc + Aq cos(2Jt), (3a) 

ujR{t) = Wc - Ao cos(2 Jt), (3b) 

where the modulation amplitude Aq is assumed to be 
much smaller than J. For simplicity, Aq and J are taken 
to be positive numbers. The sinusoidal time dependence 
of the cavity-field frequencies could be realized by modu¬ 
lating the lengths of the cavities. When the lengths of the 
left and right cavities change as L — (Ao/wc)T cos(2Jt) 
and L -|- (Ao/a;c)Tcos(2Jt), respectively, then the fre¬ 
quencies of the two cavities are modulated as in Eq. (3). 
Here L is the cavity length when the membrane and the 
outside mirrors are at rest. We assume that the motion 
of two outside mirrors is prescribed classically, and the 
back-action on such mirrors is neglected. 

To realize the modulation of the cavity lengths, one 
could use piezoelectric transducers, which convert the 
electrical signal to mechanical vibration [45]. We could 
also mount the end mirrors on mechanical resonators, 
such as cantilevers, which oscillate periodically. In ad¬ 
dition, one could shake the outside mirrors by shining 


each of them with two lasers of frequency difference 2J, 
since the beating between the two lasers yields a radia¬ 
tion pressure force modulated at a frequency 2J [46]. 

The evolution of the system state \tjj{t)) is governed by 
the Schrodinger equation H{t)\'tp{t)) = id\'tlj{t))/dt. To 
study how the frequency modulation affects the quan¬ 
tum dynamics, it is useful to introduce the transforma¬ 
tion \tp{t)) = T{t)\'4j{t)), where T{t) is a unitary operator 
defined by T{t) = Vi{t)V 2 {t)V 3 {t) with 


Vi{t) = exp[-iujct{alaL + a^aij)], 
V2{t) = expl-iJtia'laR -b a^RflL)], 


V3(t) = exp<^ -i 


- a^jOfl) -b LOmb'^b 


(4a) 

(4b) 

(4c) 


The evolution of \ipit)) is governed by the transformed 
Hamiltonian H{t) = T^{t)HT{t) — iT^{t)dT{t)/dt, which 
has the form 


m 

= ^{o^RaL - a|jai{)(e^*'^* -b e 


—4z Jt\ 


-b^[a|jai(e*('‘-^-^“)* - e-*(4J+^o)t) h.c.] 

-by (a^az, - -b 6ted<^M-e2j)t) 


-b 


% 

2 I 






(5) 


The motivation for performing the above transformation 
is that the parameters ujm and J, which are typically 
much greater than go and Aq, all appear in oscillating 
phase factors. The form of H{t) allows us to identify fast 
oscillating terms and perform the rotating-wave approx¬ 
imation. Specifically, under the conditions 

\oJm — 2J| < y, J ^ ^ 

we may retain only the slowly oscillating terms ^(a^az,— 
a|jazi)[5e-*(‘^“-2'^)* + foted<^M- 2 J)t] Eq. (5), and dis¬ 
card all other terms. This is justified within the rotating- 
wave approximation, because in order to discard those 
terms, their oscillation frequencies must be much larger 
than their corresponding coupling coefficients under con¬ 
ditions (6). Hence Eq. (5) is approximated by 

H{t) ^ y (al^Oi - aljazj)(6e-*(‘"“+ H.c.). (7) 

Such a Hamiltonian can be cast into a familiar form in 
a rotating frame defined by \tp{t)) = S(t)), with 
the unitary operator S{t) = exp{i[wc(a|^ai -b o^azz) -b 
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{u!m — 2J)b^b]t}. The transformed Hamiltonian H' = 
S^t)H{t)S{t) — iS'^{t)dS{t)/dt then reads 

H' Ri Wc(aL“i + + (wm - 2J)b% 

+ + (8) 

The right side is precisely the same form of the Hamilto¬ 
nian of the MIM optomechanical system without the tun¬ 
neling term, but with an effective mechanical frequency 
ujm — 2 J, which is significantly reduced. Consequently, 
the mechanical displacement induced by a single photon, 
which is now proportional to 5 o/[2(wm — 2J)], can be 
significant. 

As a remark, our model can be analyzed by work¬ 
ing in the Schwinger representation. Define the angu¬ 
lar momentum operators Sx = + a|jai)/2, Sy = 

i{a\f^aL — a]^aii)/2, and Sz = — a^afl)/2, then 

the Hamiltonian (2), up to a constant-of-motion term 
uJc{a\^aL + can be expressed as 

H = 2AoCOs{2Jt)Sz + 2JSx + u}Mb^b + 2gSz{b + b^), (9) 

which is the quantum Rabi model for 5=1/2 (single¬ 
photon case) and Aq = 0. The Aq term corresponds to a 
driving term. It is worth noting that the quantum Rabi 
model in the ultrastrong-coupling regime, which is diffi¬ 
cult to be realized in typical systems, can be simulated 
by reducing the effective qubit frequency [47]. 

III. SINGLE-PHOTON MECHANICAL 
DISPLACEMENT 

Let us consider an initial state of the system in which 
the membrane is in the ground state and a photon is in 
the left cavity mode, 

|V^(0)) = |^'(0)) = |l)i|0)H|0)M. (10) 

Then, by using the approximate Hamiltonian (8), we 

have 

\tjj(t)) « r(t)5(t)e"*-^'*|Vi'(0)) = |'0approx(<)) (H) 

after transforming back to the original frame. Explicitly, 
the approximate state |'(/approx(0) is given by 

IV'approx(O) = e“*®i‘ie“*‘^'=‘[cOs(Jt)|l)L|0)fl 

-*sin(Jt)|0)L|l)i?]|/3(i))M, (12) 

where |/3(t))M is a coherent state of the membrane, and 
9{t) and /3(t) are given by 

g(t] = (^- go \ ^ , go sin[(c^M - 2J)t] 

V2 4(a;M-2J)y' 4 (wm - 2J)2 ’ 

(13a) 



So^ 


FIG. 2. (Color online) (a) The fidelity F{t) versus the scaled 
time got when the parameter Ao/go takes various values, (b) 
Comparison of the mechanical displacements |(6(t))| in the 
numerical case (for Ao/go = 20, 40, and 60) with that in the 
analytical case. Other parameters are given by ujm Igo = 201 
and J/go = 100.25. 


Hence the membrane can reach a maximum displacement 
go/i^M — 2J) at the time tg = — 2J), assuming 

lom > 2J. A large displacement amplitude of the mem¬ 
brane can be obtained if 2J is close to lom- In particular, 
at 2J = the displacement amplitude in Eq. (13b) 
takes the form: j3{t) = —igotexp{—iujMt)/2, which grows 
linearly in time until it becomes too large that Hamilto¬ 
nian (2) breaks down. This is in contrast to unmodulated 
systems (i.e., Aq = 0) in which the membrane can gain 
less than one phonon for the same initial condition [48] . 

To check the validity of our approximation, we per¬ 
formed a numerical calculation of the state ['(/(t)) with 
the Hamiltonian H in Eq. (2) subjected to the frequency 
modulation defined by Eq. (3) and the same initial con¬ 
dition as above. The numerically exact ['(/(t)) is then 
compared with the ['(/approx)^)) of Eq. (12). Specifically, 
we examine the fidelity F{t) = l(V'(t)|l/’approx(t))P- 

We remark that in solving numerically the Schrddinger 
equation using the Hamiltonian (2) and the initial state 
[(/(O)) = 11 )l10)h10)m, we express the system state at 
later time as 

OO 

l^(^)) = '^[-^ni{t)\l)L\0)R +Bmi,t)\0)L\l)R]\m)A^U) 
m—0 
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Then, by the Schrbdinger equation, the probability am¬ 
plitudes Am and Bm are governed by the following set of 
coupled differential equations: 

Am{t) = - i[wc + mujM + Ao cos{2Jt)]Am{t) - iJBm{t) 
- igo[\/m+ lAm+i{t) + y/mAm-i{t)], (15a) 

Bmit) = - f[Wc -f mUJM - Ao cos{2Jt)]Bm{t) - iJAm{t) 
-I- igo[\/m,+ lBm+i{i) + y/mBm-iit)], (15b) 

which are solved numerically. The fidelity between the 
two states |'0(t)) and IV'approx(t)) can be obtained as 

Fit) = |(V'(t)|V'approx(t))P 

oo 2 

= ^^—^[coa{Jt)A*m{t)-iAn{Jt)B*m{t)\ . 

—^ \'rn\ 

m—Q ^ 

(16) 

In Fig. 2(a), we plot the fidelity as a function of the 
evolution time for various values of Ao/go- We see that a 
good fidelity [i.e., F{t) ~ 1] can be obtained for larger val¬ 
ues of Aq/^O) which is in accordance with the condition 
in Eq. (6). The fidelity curves exhibit some fast oscilla¬ 
tions, which are caused by the time factors exp(±f2Jt) 
and exp(±ia;Mt) (wmj 2J ^ qq). We also plot the time 
dependence of the mechanical displacement |(5(t))| in 
Fig. 2(b). Here the solid curve corresponds to the ana¬ 
lytical case |(6(t))| = |/3(t)|, which is independent of Aq, 
as given in Eq. (13b). The other three curves in panel 
(b) are determined by the numerical solutions with the 
corresponding values of Ap/go given in panel (a). It can 
be seen that the exact numerical results match well the 
approximate analytical result, except a slight shift of the 
period. 

Although a high ratio Ag /go would increase the fidelity 
as shown in Fig. 2, there is an upper limit of Ap accord¬ 
ing to Eq. (6). Specifically, for the approximate Hamil¬ 
tonian (7) to be valid, we also need Ap <C 16J/5. In 
other words, there is a finite range of Ap for our model 
to operate with. This is illustrated in Fig. 3 where the 
fidelity F{ts) at tg = -kIIuim — 2J) is plotted as a func¬ 
tion of lAojgo. The choice of time tg is useful because it 
is the time for the first maximum displacement. In plot¬ 
ting all the curves in Fig. 3, we use the parameters with 
wm ~ 2J = 5 p/ 2 , so that the maximum mechanical dis¬ 
placement is 2. From Fig. 3, we see that fidelity is almost 
1 in a range of Ap/gp when the conditions in Eq. (6) are 
satisfied. The range of Ap/^p giving F(ts) « 1 becomes 
wider as the ratio Jjgo increases, as illustrated by the 
blue dash-dotted curve in Fig. 3. 

IV. DISCUSSIONS 

Our analysis so far has focused on the coherent evo¬ 
lution of the modulated system. To address the effect 
of cavity field damping, we let Kc be the cavity field 
damping rate for both the left and right cavities, then 



FIG. 3. (Color online) The fidelity F{ts) at time ts versus the 
ratio Ao/go for various values of J/go, where ts = nI{ ojm—2J) 
is the time for the first peak value of |/3(t)|. For different values 
of J, the values of uim are taken to keep uim — 2J — go/2. 


with uiM — 2J ft! go/2, t ft 1/go is the time scale for the 
membrane displacement being greater than the width of 
the mechanical ground state. Therefore the condition 
go ^ i^c would ensure that the photon has a sufficient 
time to yield a significant mechanical displacement (with 
\(3\ > 1) before it leaks out of the cavity. 

More specifically, based on the approximate Hamilto¬ 
nian (8) and the initial state (10), it is not difficult to 
estimate that the phonon number (5^6) is of the order 
(;q/[4(wm — 2 J)^ -I- 4:kI] at a time long compared with the 
cavity-field lifetime, but short compared with the life¬ 
time of the mechanical oscillator, i.e., I/tm ^ t ^ 1 /kc, 
where 7 m is the damping rate of the membrane. There¬ 
fore, if go ^ Kc and com — 2J < go, then the average 
phonon number of our modulated system can be larger 
than 1, after the photon leaks out of the cavity. 

We point out that our approach can be used to cre¬ 
ate distinct mechanical superposition states [49-56] . The 
procedure consists of three steps: (i) Prepare the system 
in initial states (l/-y2)(|l)i,|0)_R ± |0)L|l)fl)|0)M- (h) 

Then let the system evolve for a suitable time, for ex¬ 
ample t = 2 tt/(u!m — 2J). (hi) Measure the state of the 
photon [44]. As long as the photon is measured in either 
the left cavity or the right cavity, then the mechanical 
resonator will be projected to a mechanical superposi¬ 
tion state. 

We note that mechanical superposed states can be 
measured by using quantum state tomography [57, 58]. 
By strongly driving another cavity mode (avoiding the 
cross talk with the single-photon state) at red sideband, 
the optomechanical coupling can be linearized into a 
beam-splitter interaction. Assume a large decay of the 
assistant cavity field such that it can be eliminated adi- 
abatically, then the state of the mechanical mode can be 
mapped to the output light [59] . As a result the mechan- 
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ical state can be known by reconstructing the state of 
the output light beam [60] . In addition, the technique of 
spectrometric reconstruction can also be used to obtain 
the state information of the mechanical mode by detect¬ 
ing the single photon spectrum [57]. 

Moreover, we would like to mention that the approach 
in this paper can also be used to enhance the optical 
nonlinearities [21]. In optomechanics, the magnitude of 
the optical nonlinearity is It can be increased by 

effectively decreasing the frequency, from wjvf to J, 

under the same coupling strength g^. 

Finally, we discuss an example that illustrates the con¬ 
ditions required for the experimental demonstration of 
the mechanism presented in this paper. The example 
consists of a vibrational mode of a silicon nitride mem¬ 
brane with an eigenfrequency ujmI‘2'^ = 1 MHz [36, 37], 
and the mechanical damping rate is less than a few hertz 
and so it may be ignored. The membrane is placed in 
the middle of a Fabry-Perot cavity of length L. A pair 
of nearly degenerate normal modes of the cavity is cho¬ 
sen such that the frequency splitting 2J Ri ojm at the 
avoided crossing. To obtain mode splitting in the mega¬ 
hertz range, the transverse modes of the cavity may be 
exploited [30, 33, 36]. Consider L = 1 mm, then go is 
about 1 kHz, which is much smaller than ujm- However, 
by applying a weak modulation with Aq = 40go ~ 40 
kHz, the effect of radiation pressure from a single pho¬ 
ton can be significantly enhanced according to our the¬ 
ory, provided that the cavity field damping rate Kc is 
smaller than gQ. Since Kc is of the order of few mega¬ 
hertz in experiments [37], this is still a main challenge to 
achieve go > Kc- Note that go > Kc is also the condition 
for the MIM optomechanical system to exhibit enhanced 
photon-photon interactions [42]. Therefore the progress 
of designs and fabrications of ultrahigh Q cavities will be 
crucial to access ultrastrong-coupling quantum effects in 
optomechanics. 

V. CONCLUSION 

In conclusion, we have proposed a method to enhance 
the mechanical effects of single photons in a two-mode 


optomechanical system by introducing a frequency mod¬ 
ulation to the cavity fields. The quantum dynamics, 
which is driven by photon hopping between the two mod¬ 
ulated cavities and optomechanical interactions, can be 
captured approximately by the Hamiltonian (8) under 
the condition (6). Such a Hamiltonian indicates that a 
single photon can displace the membrane with /3 > 1 
even though go/oJM “C 1. This method can be used to 
create distinct superposition states of the mechanical res¬ 
onator and enhance the optical nonlinearities. Finally we 
remark that these types of systems can be studied via 
quantum emulations or quantum simulations [61, 62], as 
described in [63, 64]. These describe circuit analogs of 
optomechanical systems, allowing an additional way to 
study the effects predicted here. 
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